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Abstract The availability of vector-magnctogram sequences with sufficient ac- 
curacy and cadence to estimate the temporal derivative of the magnetic field 
allows us to use Faraday's law to find an approximate solution for the electric 
field in the photosphere, using a Poloidal-Toroidal Decomposition (PTD) of the 
\Q j magnetic field and its partial time derivative. Without additional information, 

however, the electric field found from this technique is under-determined - Fara- 
day's law provides no information about the electric field that can be derived 
p/| the gradient of a scalar potential. Here, we show how additional information 



in the form of line-of-sight Dopplcr-flow measurements, and motions transverse 
to the line-of-sight determined with ad-hoc methods such as local correlation 
tracking, can be combined with the PTD solutions to provide much more ac- 
curate solutions for the solar electric field, and therefore the Poynting flux of 
electromagnetic energy in the solar photosphere. Reliable, accurate maps of the 
Poynting flux are essential for quantitative studies of the buildup of magnetic 
^ ' energy before flares and coronal mass ejections. 
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. 1. Introduction 

00 

, The launch of SDO, with its ability to measure the Sun's vector magnetic field 

anywhere on the disk with a high temporal cadence, promises to usher in a new 
era of solar astronomy. This new era of measurement demands new approaches 
for the analysis and use of this data. We show in this article how the vector 
magnetic field and Dopplcr-flow measurements that can now be made with HMI 
(jScherrer and T he HMI Team] I2005|) lead to new methods for determining the 
electric field vector, and the Poynting Flux vector 
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at the solar photosphere. The Poynting flux measures the flow of electromag- 
netic energy at the layers where the magnetic field is determined. Quantitative 
observational studies of how energy flows into the corona depend on deriving 
accurate estimates of the Poynting flux. 

Most work estimating the Sun's electric field or Poynting flux either explicitly 
or implicitly assumes that the electric field is determined by ideal MHD pro- 
cesses, and therefore the problem can be reduced to determining a velocity field 
associated with the observed magnetic-field evolution. One class of velocity esti- 
mation techniques are "Local Correlation Tracking" (LCT) methods, which es- 
sentially capture pattern motions of the line-of-sight magnetic field or white-light 
intensity. This approach was pioneered by [Novembe r and Simonl (| 1988)) . Other 
implementations include the Lockheed-Martin LCT code (jTitle et ail [1995; 
IHurlburt et aU IT995)) . "Balltracking" (|Potts, Barrett, and Diver[ 12004]) . and the 



FLCT code (jFisher and We lsch , 2008} . Another class of velocity-estimation meth- 
ods incorporate solutions of the vertical component of the magnetic induc- 
tion equation into determinations of the velocity field (jKusano et al\ [2002; 
I Welsch et al\ [20041 |Longcope[ [20041 Schuck, 2006, 2008: IChae and Sakuraill2008ll . 
The work we present in this article incorporates solutions of the three-dimensional 
magnetic induction equation, using the electric field as the fundamental variable, 
rather than the velocity field. 

The temporal evolution of the Sun's magnetic field is governed by Faraday's 
law, 

— = -V x C E . (2) 

If one can can make a map on the photosphere of dTi/dt, can one determine E by 
uncurling this equation? Addressing this question was the focus of lFisher et al.\ 
(2010}, in which a poloidal-toroidal decomposition (PTD) of the temporal deriva- 
tive of the magnetic field was used to invert Faraday's law to find E. IFisher et al.\ 
(2010) found that one could indeed find solutions for E that solve all three 
components of Faraday's law, but the solutions are not unique: the gradient of 
a scalar function can be added to the PTD solutions for E without affecting 
V x E. Fish er et al.\ (|2010p explored two different methods for determining the 
scalar function using ad-hoc and variational methods, both of which enforced the 
assumption, from ideal MHD, that E must be normal to B. Unfortunately, the 
agreement with a test case from an MHD solution, while better than conventional 
correlation-tracking methods, was still disappointing. The authors concluded 
that including additional information from other observed data was one possible 
approach for improving the electric field inversions. 

In this article, we use the same MHD simulation test case used in Welsch et al. 
(j2007[) and IFisher et aU pOTO|) to show that using Dopplcr-flow measurements 
to determine the electric scalar potential, especially in regions where the mag- 
netic field is primarily horizontal, can dramatically improve the inversion for the 
electric field and the Poynting flux. 

In Section [5] we review the PTD formalism that describes how one can derive 
the purely inductive part of the electric field from measurements that estimate 
the time derivative of B, and the technique of Section 3.2 of lFisher et a/.l(|2010[) . 
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showing how one can derive a potential electric field, which, when added to the 
inductive part of the electric field, is normal to the magnetic field. This is useful 
in generating electric-field solutions that are both consistent with Faraday's law 
and with ideal MHD, which is generally believed to be a good approximation in 
the solar photosphere. 

Section [3] argues from physical grounds why magnetic-flux emergence may 
make a large contribution to the part of the electric field attributable to a po- 
tential function. Then, starting from this argument, we derive a Poisson equation 
for an electric-field potential function that is determined primarily from knowl- 
edge of the vertical velocity field, as determined from Doppler measurements, 
and the horizontal magnetic field near polarity inversion lines where the field 
is nearly horizontal. The electric field from this contribution is then added to 
that determined from the PTD solutions. We then apply this technique to the 
MHD simulation test data, to compare the electric field from the simulation 
with that from PTD alone, and with that from combining PTD with Doppler 
measurements . 

In Sectional we try a similar approach, but instead of using contributions to 
the horizontal electric field from Doppler measurements, we use non-inductive 
contributions to the electric field determined from the FLCT correlation-tracking 
technique, applicable in regions where the magnetic field is mainly vertical. This 
technique is essentially the three-dimensional analogue of the ILCT technique 
described bv IWelsch et al.\ (|2004[) . We also try combining PTD with contribu- 
tions from both the Doppler measurements and those from FLCT, and compare 
with the simulation data. 

Our results are summarized in Section [SJ along with a discussion of where 
additional work is needed. 



2. Poloidal Toroidal Decomposition 

Here, we present only a brief synopsis of the PTD method of deriving an electric 
field E that obeys Faraday's law. More detail can be found in Section 2 of 
IFisher et all^UW\i . 

Since the three-dimensional magnetic field vector is a solenoidal quantity, one 
can express the magnetic field in terms of two scalar functions, B (the "poloidal" 
potential) and J (the "toroidal" potential), as follows: 

B = VxVx6z + Vx Ji . (3) 

Taking the partial time derivative of Equation ^ one finds 

B = V x V x Bz + V x jz . (4) 

Here, the overdot denotes a partial time derivative. We will now assume a locally 
Cartesian coordinate system, in which the directions parallel to the photosphere 
are denoted with a "horizontal" subscript h, and the vertical direction is de- 
noted with subscript z. One can then re- write Equations (J3|) and ^) in terms of 
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horizontal and vertical derivatives as 



B z , (7) 
-8.(VxB fc ), (8) 

V ft • B h . (9) 



and 



[ ^ ) + 



+ V h x jz, - W 2 h m. (6) 




One useful property of the poloidal-toroidal decomposition is that the scalar 
functions B, j, and dB/dz can all be determined by knowing the time derivative 
of the magnetic-held vector in the plane of the photosphere. By examining the 
z-component of Equation ((5J) , the z-component of the curl of Equation , and 
the horizontal divergence of Equation ©, one can derive the following three 
two-dimensional Poisson equations for B, j , and dB/dz: 

v£b = - 

V\j = -(4tt/c) X = 

and 

VlidB/dz) = 



Here, B z and denote the partial time derivatives of the vertical and horizontal 
components of the magnetic field, respectively. Solving these three Poisson equa- 
tions provides sufficient information to determine an electric field that satisfies 
Faraday's law. 

By comparing the form of Equation @ with Equations Q and © it is clear 
the following must be true: 

V x cE = -V x V x Bz - V x jz (10) 
= -V h (dB/dz)-X7 h x jz + V 2 h Bz. (11) 

Uncurling Equation ((TU)) yields this expression for the electric field E: 

cE = -V x Bz - jz - cVif) ee cE 1 - cVip. (12) 

Here, — Vtp is the contribution to the electric field from a scalar potential, for 
which solutions to Faraday's law reveal no information. The solution for E with- 
out the contribution from — V?/ 1 , E 7 , is the purely inductive solution determined 
from the PTD method. Within this article, this solution will be referred to simply 
as the PTD solution or the PTD electric field. Note that the PTD solution is not 
unique. While solutions for dB/dz are necessary to ensure that Faraday's law is 
obeyed, the PTD solution for the electric field itself depends only on B and j. 
This means that the PTD electric field is the same for distributions of B z and 
B/j which have differing values of • B^, but the same values of (V^ x B^,) ■ z 
and B z . Thus the PTD solutions for E 1 arc under-determined. 
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IFisher et al.\ l|2010p described two techniques for deriving an electric- field con- 
tribution from a scalar potential, in an effort to resolve the under-determined 
nature of the PTD solutions. The first technique, described in Section 3.2 of that 
article, presents an ad-hoc iterative method for deriving a scalar potential electric 
field which, when added to the PTD solution, results in an electric field that is 
normal to B, and hence consistent with ideal MHD. The second technique, based 
on a variational method, finds a scalar potential electric field that, when added 
to the PTD solution, minimizes the area integral of |E| 2 or |v| 2 . When compared 
to the original electric field from the simulation test case, the iterative method 
applied to the PTD solutions showed a qualitative consistency, but not detailed 
agreement with the simulation electric field, while the electric field computed 
with the variational technique showed poor agreement. IFisher et aZ.I (|2010f) con- 
cluded that significant improvement in the agreement of the inverted electric field 
with the real electric field requires additional observational information beyond 
the temporal evolution of B. 



3. The Importance of Doppler Flow Measurements to the Electric 
Field 

We argue here that when flux emergence occurs, much of the missing infor- 
mation about non-inductive contributions to the electric field is contained in 
Doppler flow information (see also Ravindra, Longcope, and Abbett| 120081 and 



Sch uck et aU I2010|) . particularly near polarity inversion lines (PILs), where the 
horizontal magnetic field is much stronger than the vertical field. We illustrate 
this point with a simple thought-experiment, shown schematically in Figure [TJ 
Consider the emergence of new magnetic flux in an idealized bipolar-flux system, 
where the PIL maintains its orientation as flux continuously emerges from below 
the photosphere. Imagine that vector magnctogram and Doppler observations are 
taken from a vantage point normal to the solar surface. Let us focus attention 
on what is happening near the center of the PIL. Suppose the magnetic field 
there remains time-invariant as flux continues to emerge, so the time derivative 
of the magnetic field there is zero, implying that Faraday's law cannot be used 
to infer the physics of the emerging flux. Yet the electric field at this location 
should be very large, driven by the upward motion of the plasma carrying the 
strong, horizontal field. In this case, magnetic-flux emergence will have a strong 
inductive signature at the edges of the idealized active region, where the vertical 
magnetic field is changing rapidly, but not near the center of the PIL. Thus, it 
seems plausible that the electric field near PILs in more realistic emerging-flux 
configurations will have a significant non-inductive component. 

Starting from this perspective, we have explored enhancements to the PTD 
method that use Doppler flow information to more tightly constrain the PTD 
electric field solutions, with the additional assumption that the photospheric 
electric field is primarily governed by ideal MHD processes. Directly above PILs, 
the vertical velocity and the observed horizontal component of the magnetic field 
unambiguously determine the horizontal electric field: 

cEf = -v z z x B h , (13) 
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Figure 1. Schematic illustration of the emergence of new flux over a time interval At, viewed 
in a vertical plane normal to the polarity inversion line (PIL) in an idealized bipolar flux system. 
The emerging flux is rising at a speed v z , which could be inferred by the Dopplcr shift measured 
by an observer viewing the PIL from above. The width of the bipolar flux system (the distance 
from the outer edge of one pole to the outer edge of the other pole) at the beginning of At is 
2xo- Notice that the change in B z at the outer edges of the emerging flux region is large, while 
the change in B z at the PIL itself — where the flux is actually emerging — is zero (see text). 
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where we assume that |B Z |/|S^| is small. If we can use line-of-sight Doppler 
velocity measurements to estimate v z , we add a powerful constraint to the 
PTD solution for the electric field. Of course, we would like to use the Doppler 
information away from PILs as well, but are hindered by two complications: 
i) flows parallel to the magnetic field will not affect the electric field at all, 
but may contribute to the observed Doppler velocity signal, and it) when the 
vertical component of the magnetic field [B z ] becomes significant compared to 
the horizontal field [Bh], there is an additional contribution to the horizontal 
electric field from flow parallel to the surface, which is not accounted for. 

We now develop a formal solution for a non-inductive contribution to the 
electric field that includes information from Doppler-shift measurements, and 
apply it to a test case with a known electric field. First, from the pair of synthetic 
vector magnetograms taken from the ANMHD simulation test case described in 
IWelsch et all ([2007)) and Section 3.1 of IFisher et ail QZGTty . we use the PTD 
method to find an electric-field solution, neglecting any contribution from a 
scalar electric-field potential function. We use the numerical techniques and 
boundary conditions described in Section 3.1 of lFisher et aZ.I(|2010p . Second, we 
compute a candidate horizontal electric field from vertical velocities taken from 
the simulation as synthetic Doppler-flow measurements, and horizontal magnetic 
fields from the synthetic vector magnetograms, from Equation (|13[) above. This 
electric field is then multiplied by a "confidence function" , which is near unity 
at PILs, but decreases to zero when \B Z /Bh\ is no longer small. This reflects 
our lack of confidence in the accuracy of this horizontal Doppler electric field 
in those locations, for the reasons described earlier. The specific form for the 
confidence function is probably not important. Here, we assume the confidence 
function w is given by 

w = cxp[-(\B z \/\B h \) 2 /a 2 ] , (14) 

where a is a free parameter that can be adjusted, and in the specific cases shown 
in this article was set somewhat arbitrarily to 0.6. We define the "modulated" 
electric field within the plane of the magnetogram as 

Ef=«< (15) 

Third, we take the divergence of this modulated horizontal electric field E^ 7 , 
and find the electric-potential function that can best represent it by setting 

cE* = -V hX , (16) 

where x solves the Poisson equation 

V£x = -V h • cEf . (17) 

Because the synthetic vector magnetograms and Doppler flows taken from the 
MHD simulation use periodic boundary conditions, we use FFT techniques to 
solve Equation ([T7|) . Adding this contribution onto the PTD solutions means 
that information about the electric field at PILs has been incorporated, while 
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also maintaining consistency with Faraday's Law, since E x has no curl. Since 
we generally expect ideal MHD to be a good approximation for conditions in 
the solar photosphere, we then remove the components of E parallel to B by 
adding the electric field from a second potential function ip, using the iterative 
technique described in Section 3.2 of lFisher et al.\ (|2010p : 

cE tot = cE 7 + cE x - VV; , (18) 

where V?/> • B = (cE 1 + cE x ) • B. 

The resulting solutions for E are shown in the third row of Figure [51 with 
a scatterplot comparison of S z of the PTD method and the PTD plus Doppler 
information with the actual simulation electric fields shown in the top two panels 
of Figure [31 These portions of the figures show that the recovery of the electric- 
field components and the Poynting flux is dramatically better than PTD alone. 

4. How Important are Horizontal, Non-inductive Flows? 

In the previous section, we considered the role of Doppler-flow measurements 
in determining non-inductive contributions to the horizontal electric field, and 
found that combining this information with the PTD solutions for Faraday's law 
results in a dramatic improvement in the recovery of the electric field. However, 
this treatment neglects possible contributions to the horizontal electric field away 
from PILs where a cross product of horizontal velocity with vertical magnetic 
field could also contribute to the horizontal electric field. Contributions to the 
horizontal electric field that solve the induction equation have already been 
incorporated by the PTD solutions, but as with vertical velocities, there could 
be a sub-space of horizontal flows that do not contribute to Faraday's law. 

To evaluate this effect, we estimate horizontal velocities using the FLCT local- 
correlation tracking (LCT) code (see IFisher and Welschl I2008P , available from 
|http: / /solarm uri.s sl.berkeley.edu/~fisher/ public/software/FLCT/C_VERSIONS/ us- 
ing images of the vertical component of the magnetic field. Velocities were not 
computed for pixels with a vertical magnetic field strength below 370G (see 
discussion in Wc lsch et all I2007P , with the windowing parameter a set to five 
pixels. The low-pass filtering option was not invoked. The result is a map of the 
apparent horizontal-velocity field [U^ = U x x + U y y] computed at the strong 
vertical magnetic-field locations, and with velocities at all other locations set to 
zero. A candidate horizontal electric field is estimated by setting 

cE lct = _ Vh x ~ Bz ( 19 -j 

To consider only non-inductive contributions from TJh, we perform the same 
general operation as in the previous section, namely to multiply E^ CT by a con- 
fidence function, and then eliminate the inductive part of the electric field. Here, 
the confidence function will be the complement of the confidence function used 
for the Doppler case, since the LCT estimates are nearly useless near PILs, where 
the Doppler results should be reliable, while the LCT results should be best when 
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Figure 2. Top row: The three components of the electric field and the vertical Poynting flux 
from the MHD reference simulation of emerging magnetic flux in a turbulent convection zone. 
Second row: The inductive components of E and ,S Z determined using the PTD method. Third 
row: E and S z derived by incorporating Dopplcr flows around PILs into the PTD solutions. 
Note the dramatic improvement in the estimate of S z . Fourth row: E and S z derived by 
incorporating only non-inductive FLCT derived flows into the PTD solutions. Note the poorer 
recovery of E x , E y , and S z relative to the case that included only Dopplcr flows. Fifth row: 
E and S z derived by including both Doppler flows and non-inductive FLCT flows into the 
PTD solutions. Note the good recovery of E x , E y , and S z , and the reduction in artifacts in 
the low-field regions for E y (best viewed in the electronic version of the article). 
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Figure 3. Upper left: A comparison of the vertical component of the Poynting flux derived 
from the PTD method alone with the actual Poynting flux of the MHD reference simulation. 
Upper right: A comparison between the simulated results and the improved technique that 
incorporates information about the vertical flow field around PILs into the PTD solutions. 
Lower left: Comparison of the vertical Poynting flux when non-inductive FLCT-derived flows 
are incorporated into the PTD solutions. Lower right: Comparison of the vertical Poynting flux 
when both Doppler flow information and non-inductive FLCT-derived flows are incorporated 
into the PTD solutions. Each scatterplot also shows the computed linear correlation coefficient, 
as well as the slope of the fit derived with IDL's LADFIT function. Poynting flux units are in 
[10 5 G 2 km s- 1 ] 



the magnetic field is mostly vertical (and where the Doppler measurements are 
useless). 

We define cE^ = — V/i£, and assume that 



(1 - iu)c^ 



LCT 



(20) 



Once this equation has been solved and E^ has been computed, it can be 
added to the PTD solutions for E 7 , and as in the previous section, a second 
potential solution can be found that eliminates components of E parallel to 
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B. Note that combining the PTD solutions with in this way is like the 
approach used in the ILCT technique described by IWelsch et al. \ (|2004p , except 
that solutions of a single component of the induction equation are replaced by 
solutions to all three components of the induction equation. 

The resulting electric field and Poynting flux can be compared to the actual 
case, the un-altcred PTD case, and the case where only the Dopplcr information 
is used. The electric field and Poynting flux results are shown as the fourth row 
of panels in Figure [2] and a scatterplot of the Poynting flux values with the actual 
values is shown in the lower left panel of Figure [3J While the overall performance 
of the FLCT case is better than that of PTD alone, it is not significantly better 
than simply applying the iterative method directly to the PTD results as was 
described in Section 3.2 of IFisher et al.\ (|2010[) . It is definitely not as good as 
the performance we show from the Doppler-only case. We conclude that most 
of the useful information about the non-inductive electric field, at least for this 
particular simulation of strong flux emergence, is contained within the Doppler 
flow information. 

Does the LCT information, when added to the Doppler-flow information, 
significantly improve the resulting estimate for the electric field? To answer this 
question, we have added both the LCT and Doppler electric-field information to 
the PTD solutions, and again found a potential function to eliminate components 
of E parallel to B. The resulting electric field and Poynting- flux maps are shown 
in the fifth row of panels in Figure [2 and a scatterplot of the vertical Poynting 
flux is shown in the lower-right panel of Figure [3J 

The linear correlation coefficient in the Poynting-flux scatterplot is not signfi- 
cantly improved by adding the LCT results to the Doppler results, but the slope 
of the fit (determined by using IDL's LADFIT function) is somewhat better. 
Further, examining the maps of E x and E y show a reduction in artifacts in the 
behavior of the recovered electric-field components, compared to the Doppler 
and LCT cases. We conclude that at least for this simulation, which exhibited 
strong flux emergence, most of the additional useful informat ion beyond solutions 
to Faraday's law is contained within the Doppler velocity measurements, with 
some additional improvement when non-inductive LCT-dcrivcd electric fields arc 
added. 

Finally, we wish to add a comment about solutions to the PTD equations 
themselves. The PTD solutions used in this article did not use FFT solutions for 
B and J , even though the simulations are periodic, but instead used Neumann 
boundary conditions for B for the reasons described in Section 2.2 of IFisher et al.\ 
(2010]). For the current study, we compared the results of using FFT solutions 
of the PTD equations with those shown in the figures in this article, and found 
noticeable degradations in the fits of the model Poynting fluxes to the actual 
model values. If one is interested in the most accurate reconstruction of the 
vertical Poynting flux, we recommend not using FFT solutions of the PTD 
equations. 
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5. Discussion and Conclusions 

We have reviewed how the PTD solutions of Faraday's law for E can be found 
using temporal sequences of vector magnetograms that can be obtained with the 
HMI instrument on NASA's SDO mission. We discussed why these solutions are 
under-determined, and the importance of determining the contributions to the 
electric field that can be derived from a scalar potential. 

We demonstrate, using simulation data where the true electric field is known, 
that knowledge of the vertical-velocity field (obtainable by Doppler measure- 
ments) can provide important information about the electric field. When this 
information is combined with the PTD solutions of Faraday's law, dramatically 
more accurate recovery of the true electric field is possible. We find that ad- 
ditional information about flows from local correlation-tracking methods can 
also be combined with the PTD solutions, but the additional information is 
signficantly less important than that from the Doppler measurements. We are 
able to quantitatively reconstruct the electromagnetic Poynting flux in the sim- 
ulations by using our combination of the PTD solutions and those from Doppler 
measurements. 

This "proof-of-concept" demonstration argues strongly for the development of 
electric-field and Poynting-flux tools to be used routinely in the analysis of HMI 
vector magnetic-field measurements. Routinely available Poynting-flux maps will 
be useful for scientific studies of flare-energy buildup, understanding the flow of 
magnetic energy in the solar atmosphere prior to CME initiation, and will aid 
in understanding the flow of energy that heats the corona. Further, the PTD 
formalism for the magnetic field itself (Equation ((SJ) allows for a straight-forward 
decomposition of the Poynting flux into changes in the potential-field energy, and 
the flux of free magnetic energy (see I Wclsch , 20011 and the end of Section 2.1 of 
Fish er et al.\ I2010p . The flux of free magnetic energy is especially important in 
determining how energy builds up in flare-productive active regions. 

To find solutions for E and the Poynting flux S using the PTD formalism 
plus Doppler measurements requires only the solution of three two-dimensional 
Poisson equations. While real vector-magnetogram patches will not have periodic 
boundary conditions (as were employed in this article), straightforward numeri- 
cal techniques exist to solve these equations routinely. Preliminary investigations 
also indicate that generalizing the PTD solutions and Doppler measurements 
to cases of non-normal viewing angle will be straightforward. In our opinion, 
the major obstacle that remains before such solutions can be routinely applied 
to the HMI data, is a detailed understanding of how measurement errors and 
disambiguation errors in the vector magnetograms will affect the solutions, and 
how the effects of these errors arc best ameliorated. 

Acknowledgements This research was funded by the NASA Hcliophysics Theory Program 
(grant NNX08AI56G), the NASA Living-With-a-Star TR&T Program (grant NNX08AQ30G), 
by the NSF SHINE program (grants ATM0551084 and ATM0752597), and support from NSF's 
AGS Program (grant ATM0641303) for our participation in the University of Michigan's 
CCHM Project. The authors are grateful to US taxpayers for providing the funds necessary to 
perform this work. The authors wish to acknowledge Dick Canfield for his pioneering work in 



SOLA: ms_20110315_sp.tex; 18 March 2011; 1:06; p. 12 



Electric Fields and Poynting Fluxes 



the use of vector magnetograms in solar physics. The inspiration for the work described here 
can be traced to a Solar MURI workshop held at UC Berkeley in 2002, in which Dick Canficld 
played a major role in defining long-term research goals for the use of vector magnetograms 
in quantitative models of the Sun's atmosphere. 



References 

Chae, J., Sakurai, T.: 2008, A Test of Three Optical Fl ow Techniques-LCT, DAVE, and NAVE. 

Astrophys. J. 689, 593-612. doi: 10. 1086/592761] 
Fisher, G.H., Welsch, B.T.: 2008, FLCT: A Fast, Efficient Method for Performing Local Cor- 
relation Tracking. In: Howe, R., Komm, R.W., Balasubramaniam, K.S., Petrie, G.J.D. 

(eds.) Subsurface and Atmospheric Influences on Solar Activity, CS-383, Astronom. Soc. 

Pacific, San Francisco, 373-380 (also arXiv:0712.4289). 
Fisher, G.H., Welsch, B.T., Abbett, W.P., Bercik, D.J.: 2010, Estimating Elec- 
tric Fields from Vector Magnctogram Sequences. Astrophys. J. 715, 242-259. 

doi:10.1088/0004-637X/715/l/242 
Hurlburt, N.E., Schrijver, C.J., Shine, R.A., Title, A.M.: 1995, Simulated MDI Observations of 

Convection. In: Hoeksema, J.T., Domingo, V., Fleck, B., Battrick, B. (eds.) Helio seismology 

SP-376, ESA, Noordwijk, 239. 
Kusano, K., Maeshiro, T., Yokoyama, T., Sakurai, T.: 2002, Measurement of magnetic hclicity 

injection and free energy loading into the solar corona. Astrophys. J. 577, 501-512. 
Longcope, D.W.: 2004, Inferring a photosphcric velocity field from a sequence of vector 

magnetograms: the minimum energy fit. Astrophys. J. 612, 1181-1192. 
November, L.J., Simon, G.W.: 1988, Precise proper-motion measurement of solar granulation. 

Astrophys. J. 333, 427-442. 
Potts, H.E., Barrett, R.K., Diver, D.A.: 2004, Balltracking: An highly effici ent method fo r 

tracking flow fields. Astron. Astrophys. 424, 253-262. doi 10.1051/0004-6361:20035891 
Ravindra, B., Longcope, D.W., Abbett, W.P.: 2008, Inferring Photosphcric Velocity Fields 

Using a Combination of Minimum Energy Fit, Loca l Correlation Tracking, and Doppler 

Velocity. Astrophys. J. 677, 751-768. doi: 10.1086/528363 
Schcrrer, P.H., The HMI Team: 2005, The helioscismic and magnetic imager for the solar 

dynamics observatory. AGU Spring Meeting Abstracts, 43-05. 
Schuck, P.W.: 2006, Tracking Magnetic Footpoints with the Magnetic Induction Equation. 

Astrophys. J. 646, 1358-1391. doi 10.1086/505015 
Schuck, P.W.: 2008, Tracking Vector Magnetograms with the Magnetic Induction Equation. 

Astrophys. J. 683, 1134-1152. doi 10.1086/589434 
Schuck, P.W., Sun, X., Muglach, K., Hoeksema, J.T.: 2010, Tracking Vector Magnetograms 

from the Solar Dynamics Observatory. AGU Fall Meeting Abstracts, A7. 
Title, A.M., Hurlburt, N.E., Schrijver, C.J., Shine, R.A., Tarbcll, T.: 1995, Observations of 

Convection. In: Hoeksema, J.T., Domingo, V., Fleck, B., Battrick, B. (eds.) Helioseismology 

SP-376, ESA, Noordwijk, 113. 
Welsch, B.T.: 2006, Magnetic Flux Cancellation and Coronal Magnetic Energy. Astrophys. J. 

638, 1101-1109. doi 10.1086/498638 
Welsch, B.T., Fisher, G.H., Abbett, W.P., Regnier, S.: 2004, ILCT: Recovering Photo- 
sphcric Velocities from Magnetograms by Combining the Induction Equation with Local 

Correlation Tracking. Astrophys. J. 610, 1148-1156. doi 10.1086/421767] 
Welsch, B.T., Abbett, W.P., DeRosa, M.L., Fisher, G.H., Georgoulis, K., Kusano, M.K., 

Longcope, D.W., Ravindra, B., Schuck, P.W.: 2007, Tests and comparisons of velocity 

inversion techniques. Astrophys. J. 670, 1434-1452. 



SOLA: ms_20110315_sp.tex; 18 March 2011; 1:06; p. 13 



SOLA: ms_20110315_sp.tex; 18 March 2011; 1:06; p. 14 



